The role of post-failure brittleness of soft rocks in the assessment of stability of intact masses: FDEM technique applications to ideal problems

Strain-softening under low confinement stress, i.e. the drop of strength that occurs in the post-failure stage, represents a key factor of the stress-strain behavior of rocks. However, this feature of the rock behavior is generally underestimated or even neglected in the assessment of boundary value problems of intact soft rock masses. This is typically the case when the stability of intact rock masses is treated by means of limit equilibrium or finite element analyses, for which rigid-plastic or elastic perfectly-plastic constitutive models, generally implementing peak strength conditions of the rock, are respectively used. In fact, the aforementioned numerical techniques are characterized by intrinsic limitations that do not allow to account for material brittleness, either for the method assumptions or due to numerical stability problems, as for the case of the finite element method, unless sophisticated regularization techniques are implemented. However, for those problems that concern the stability of intact soft rock masses at low stress levels, as for example the stability of shallow underground caves or that of rock slopes, the brittle stress-strain response of rock in the post-failure stage cannot be disregarded due to the risk of overestimation of the stability factor. This work is aimed at highlighting the role of post-peak brittleness of soft rocks in the analysis of specific ideal problems by means of the use of a hybrid finite-discrete element technique (FDEM) that allows for the simulation of the rock stress-strain brittle behavior in a proper way. In particular, the stability of two ideal cases, represented by a shallow underground rectangular cave and a vertical cliff, has been analyzed by implementing a post-peak brittle behavior of the rock and the comparison with a non-brittle response of the rock mass is also explored. To this purpose, the mechanical behavior of a soft calcarenite belonging to the Calcarenite di Gravina formation, extensively outcropping in Puglia (Southern Italy), and the corresponding features of the post-peak behavior as measured in the laboratory, have been used as a reference in this work, as well as the typical geometrical features of underground cavities and rock cliffs, as observed in Southern Italy, have been adopted for the simulations. The numerical results indicate the strong impact for the assessment of stability when rock post-peak brittleness is accounted for, if compared with perfectly plastic assumptions, and the need for adopting numerical techniques, as the FDEM approach, to take properly into account this important aspect of the rock behavior is highlighted

A Geomechanical Approach to Landslide Hazard Assessment: The Multiscalar Method for Landslide Mitigation

AbstractThe landslide hazard assessment, when based on the deterministic diagnosis of the processes, can be pursued only through the interpretation and the geo-hydro-mechanical modelling of the slope equilibrium. In practice, though, landslide hazard assessment is still seldom dealt with slope modelling, in particular when it addresses vast areas, where either heuristic or statistical methods do not entail any geo-hydro-mechanical knowledge of slope features and stability. The Multiscalar Method for Landslide Mitigation (MMLM) is an original methodological approach for intermediate to regional landslide hazard assessment. It is based on the geo-hydro-mechanical knowledge achieved from the application of a stage-wise diagnostic methodology of the landslide mechanism at the slope scale. The paper discusses the main steps of the MMLM aiming at diagnoses of landslide hazard based on hydro-mechanics, for small scale hazard mapping (at the large area)

Assessment of post-failure evolution of a large earthflow through field monitoring and numerical modelling

AbstractThe analysis of the residual hazard existing after the emergency phases generated by the activation or reactivation of landslides is rarely taken into account in a proper manner. However, the assessment of landslide post-failure evolution should represent a key factor to control potential landslide reactivations and prevent new landslide-induced damages. This paper presents the results of a long-term field monitoring activity performed in the years after the emergency phase of the Montaguto (Italy) earthflow reactivation occurred in 2010 as well as the results of 2-D and 3-D numerical analyses aimed at interpreting the post-emergency landslide behaviour. The results of the numerical simulations, which agree well with the in situ monitoring data, allow to define a conceptual model of the earthflow behaviour that is related to the pore water pressure variations resulting from the drained or undrained processes occurring in the landslide body. The study proposed confirms a general reduction of the landslide activity, as well as allows to detect the factors that control the residual activity existing in a specific area of the landslide and to infer possible critical scenarios for landslide reactivations

An Integrated Procedure to Assess the Stability of Coastal Rocky Cliffs: From UAV Close-Range Photogrammetry to Geomechanical Finite Element Modeling

The present paper explores the combination of unmanned aerial vehicle (UAV) photogrammetry and three-dimensional geomechanical modeling in the investigation of instability processes of long sectors of coastal rocky cliffs. The need of a reliable and detailed reconstruction of the geometry of the cliff surfaces, beside the geomechanical characterization of the rock materials, could represent a very challenging requirement for sub-vertical coastal cliffs overlooking the sea. Very often, no information could be acquired by alternative surveying methodologies, due to the absence of vantage points, and the fieldwork could pose a risk for personnel. The case study is represented by a 600 m long sea cliff located at Sant\u2019Andrea (Melendugno, Apulia, Italy). The cliff is characterized by a very complex geometrical setting, with a suggestive alternation of 10 to 20 m high vertical walls, with frequent caves, arches and rock-stacks. Initially, the rocky cliff surface was reconstructed at very fine spatial resolution from the combination of nadir and oblique images acquired by unmanned aerial vehicles. Successively, a limited area has been selected for further investigation. In particular, data refinement/decimation procedure has been assessed to find a convenient three-dimensional model to be used in the finite element geomechanical modeling without loss of information on the surface complexity. Finally, to test integrated procedure, the potential modes of failure of such sector of the investigated cliff were achieved. Results indicate that the most likely failure mechanism along the sea cliff examined is represented by the possible propagation of shear fractures or tensile failures along concave cliff portions or over-hanging due to previous collapses or erosion of the underlying rock volumes. The proposed approach to the investigation of coastal cliff stability has proven to be a possible and flexible tool in the rapid and highly-automated investigation of hazards to slope failure in coastal areas

The neglected role of karst features in rock mass characterization and stability assessment

Stability analyses in karst settings, whether to assess the equilibrium conditions of natural slopes or to design engineering interventions, coexist with a significant uncertainty related to difficulties in modelling karst features. As a matter of fact, most of the rock mass classification systems do not directly take into account the presence of karst structures such as voids, conduits or caves, which can strongly influence the mechanical behaviour and the water flow in rock masses. In the last decades, the identification and characterization of discontinuity systems for rock mass characterization, aimed at stability analyses, have been intensively investigated by means of remote sensing techniques. However, semi-automatic or automatic methods for the extraction of discontinuities from point clouds are not easily applicable in karst because surface and subsurface features produce irregular surfaces, which are difficult to classify even using the most-advanced algorithms. This occurs even more heavily in the case of soft rocks, such as calcarenites. In this study, a demonstration of the influence of karst features in rock mass characterization and slope stability assessment is presented. First, the results of the Discontinuity Set Extractor (DSE) software used on an appropriate case study show that the irregular surfaces produced by carbonate dissolution, further enhanced by weathering, caused an incorrect classification of the discontinuity sets. Second, a high-resolution Digital Outcrop Model (DOM) was used to generate a very fine mesh (average element size = 35 cm, to take into account the large-scale karst structures) and to carry out 3-D numerical stability analyses by means of Finite Element Method, using a continuum-based approach. Although in the current conditions the examined slope is stable, the results illustrate that the maximum shear strain is localized in correspondence of the karst features (e.g. caves and voids) and at the sea level. By applying the Shear Strength Reduction method, it was found out that weathering processes can cause the same structures to be under yield and lead to localized failures. In addition, the key role that the discontinuities (extracted using an ad-hoc procedure) play on the rock mass mechanical behaviour was investigated using a 2-D FEM, based on a discontinuum approach. The results, which are in agreement with field observations, point out that karst processes, which features are characterized by the highest values in pervasiveness and aperture of the discontinuity systems and tend to reduce the rock bridges over time, need to be implemented in the rock mass classification systems and in numerical modelling techniques to avoid incorrect results

A multidisciplinary investigation of deep-seated landslide reactivation triggered by an extreme rainfall event: a case study of the Monesi di Mendatica landslide, Ligurian Alps

AbstractIn November 2016, an extreme rainfall event affected the Ligurian Alps (NW Italy). Consequently, several landslides and debris flows occurred in the upper Tanarello stream basin. In particular, the village of Monesi di Mendatica was severely damaged by two landslide phenomena: the activation of a rotational landslide, which caused the total collapse of two buildings and part of the main road, and the reactivation of a deep-seated planar massive and a complex landslide, which widely fractured most of the buildings in the village. The latter phenomenon was mostly unknown and had never been monitored prior to the 2016 event. Due to the extensive damage, the village of Monesi was completely evacuated, and the road connecting a ski resort area in the upper part of the valley was closed. Furthermore, a potentially dangerous situation related to the eventual progressive evolution of this landslide that could cause a temporary occlusion of the Tanarello stream still remains. For this reason, we defined the landslide behaviour, triggering conditions and chronological evolution leading to the 2016 event using a multidisciplinary approach. This approach consisted of field surveys, satellite DInSAR time series analyses, digital image correlation techniques, rainfall records analyses, postevent monitoring campaigns and subsurface investigation data analyses, and numerical modelling. This multidisciplinary approach enhanced our understanding of this landslide, which is fundamental to better comprehend its behaviour and possible evolution

Multi-scale approach to analyse the evolution of soft rock coastal cliffs and role of controlling factors: a case study in South-Eastern Italy

The evolution of soft rock coastlines is strictly related to natural and anthropogenic conditions, which in some cases can determine also an acceleration of coastal retreat. The recent evolution of a soft rocky coastal stretch on the Italian southern Adriatic sector is analysed. To investigate the most important contributing factors to coastal evolution, a detailed multi-temporal (1954–2017) morphological analysis with photo-interpretation is carried out, which allows gaining a general understanding of the cliff evolution. The coastal retreat was then evaluated by exploiting the Digital Shoreline Analysis System (DSAS) tool. From that, site-specific analyses are carried out to explore the role of four environmental factors: rock mechanical properties, geostructural setup of the rock mass along the coast, emerged rock platform at the toe, and the storm waves. Analyses highlighted how the most significant retreat is detected along cliff segments characterized by the presence of low strength rocks susceptible to water-induced weakening, moderate fracturing degree of the rock mass, absence of the emerged platform at the toe, and wave storms. This work shows how a multi-scale methodology could represent an efficient approach to gain an interpretation of the instability processes and thus prepare risk mitigation plans and land management strategies

Numerical modelling of slope–vegetation–atmosphere interaction: an overview

The behaviour of natural and artificial slopes is controlled by their thermo-hydro-mechanical conditions and by soil–vegetation–atmosphere interaction. Porewater pressure changes within a slope related to variable meteorological settings have been shown to be able to induce soil erosion, shrinkage–swelling and cracking, thus leading to an overall decrease of the available soil strength with depth and, ultimately, to a progressive slope collapse. In terms of numerical modelling, the stability analysis of partially saturated slopes is a complex problem and a wide range of approaches from simple limit equilibrium solutions to advanced numerical analyses have been proposed in the literature. The more advanced approaches, although more rigorous, require input data such as the soil water retention curve and the hydraulic conductivity function, which are difficult to obtain in some cases. The quantification of the effects of future climate scenarios represents an additional challenge in forecasting slope–atmosphere interaction processes. This paper presents a review of real and ideal case histories regarding the numerical analysis of natural and artificial slopes subjected to different types of climatic perturbations. The limits and benefits of the different numerical approaches adopted are discussed and some general modelling recommendations are addressed